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We quantitatively analyze the dynamics of the quantum phase distribution associated with the 
reduced density matrix of a system, as the system evolves under the influence of its environment 
with an energy-preserving quantum nondemolition (QND) type of coupling. We take the system to 
be either an oscillator (harmonic or anharmonic) or a two- level atom (or equivalently, a spin- 1/2 
system), and model the environment as a bath of harmonic oscillators, initially in a general squeezed 
C"*** ■ thermal state. The impact of the different environmental parameters is explicitly brought out as 

the system starts out in various initial states. The results are applicable to a variety of physical 
systems now studied experimentally with QND measurements. 

£H ■ PACS numbers: 03.65.Yz, 03.65.Vf, 42.50.Ct 
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I. INTRODUCTION 



The theory of open quantum systems addresses the problems of damping and dephasing in quantum systems by its 
■ assertion that all real systems of interest are in fact 'open' systems, each surrounded by its environment. Quantum 
optics provided one of the first testing grounds for the application of the formalism of open quantum systems 1|. 
Application of open system ideas to other areas of physics was intensified by the works of Caldeira and Leggett 2 ] , 
and Zurek Q among others. Most such studies are based on a model describing quantum Brownian motion of a simple 
harmonic oscillator in a harmonic oscillator environment. In such a model studied by Caldeira and Leggett @, the 
coordinate of the particle was coupled linearly to the harmonic oscillator reservoir, and it was also assumed that the 
system and the environment were initially separable. The treatment of the quantum Brownian motion has since been 
generalized to the physically reasonable initial condition of a mixed state of the system and its environment by Hakim 
and Ambegaokar [J], Smith and Caldeira (5], Grabert, Schramm and Ingold [6[, and by us for the case of a system in a 
Stern-Gerlach potential Q , and also for the quantum Brownian motion with nonlinear system-environment couplings 

The recent upsurge of interest in the problem of open quantum systems is because of the spectacular progress in 
manipulation of quantum states of matter (atoms, or bosonic or fermionic gases or molecules), encoding, transmission 
and processing of quantum information, for all of which understanding and control of the environmental impact are 
essential. For such open quantum systems, there exists an important class of energy-preserving measurements in which 
dephasing occurs without damping of the system. This may be achieved with a particular type of coupling between 
the system and its environment, viz., when the Hamiltonian H$ of the system commutes with the Hamiltonian Hsr 
describing the system-reservoir interaction, i.e., Hsr is a constant of motion generated by Hs 0, E3, [Hj|. This 
condition describes a particular type of quantum nondemolition (QND) measurement scheme. 

In general, a class of observables that may be measured repeatedly with arbitrary precision, with the influence 
of the measurement apparatus on the system being confined strictly to the conjugate observables, is called QND or 
back-action evasive observables 1 1 2L 1 1 3B . Such a measurement scheme was originally introduced in the context of the 
detection of gravitational waves [l4, Il5l. ITfil. Il7l. Il8l. Il9l. l2fjj . The experimental progress on QND measurements has been 
summarized in the review [2l| and the dynamics of decoherence in continuous atom-optical QND measurements studied 
in [22l ]. In addition to its relevance in ultrasensitive measurements, a QND scheme provides a way to prepare quantum 
mechanical states which may otherwise be difficult to create, such as Fock states with a specific number of particles. 
It has been shown that the accuracy of atomic interferometry can be improved by using QND measurements of the 
atomic populations at the inputs to the interferometer [23| . QND systems have also been proposed for engineering 
quantum dynamical evolution of a system with the help of a quantum meter [24[. We have recently studied such 
QND open system Hamiltonians for two different models of the environment describable as baths of either oscillators 
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or spins, and found an interesting connection between the energy-preserving QND Hamiltonians and the phase space 
area- preserving canonical transformations [25j . 

As stated above, in the context of energy-preserving QND systems, the only effect of the environment on the system 
is dephasing and it is a natural question to ask about the pattern of diffusion of 'phases' in such a situation. Such a 
question is particularly relevant in the context of a number of practical phase measurement schemes I23l. |26l . 

What is t he p recise meaning of the quantum mechanical phase? The quantum description of phases [27[ has a long 
history [28|, |29| . l30l . l3ll |32| . Pegg and Barnett [3l[ , following Dirac [28| , carried out a polar decomposition of the 
annihilation operator and defined a hermitian phase operator in a finite-dimensional Hilbert space. In their scheme, 
the expectation value of a function of the phase operator is first carried out in a finite-dimensional Hilbert space, and 
then the dimension is taken to the limit of infinity. However, it is not possible to interpret this expectation value as 
that of a function of a hermitian phase operator in an infinite-dimensional Hilbert space 1331 1341 ■ To circumvent this 
problem, the concept of phase distribution for the quantum phase has been introduced [33l. l3a|. In this scheme, one 
associates a phase distribution to a given state such that the average of a function of the phase operator in the state, 
computed with the phase distribution, reproduces the results of Pegg and Barnett. 

In this paper we address the problem of quantum phase diffusion and study the dynamics of the quantum phase dis- 
tribution associated with the reduced density matrix of the system for a number of situations of practical importance, 
as the system evolves under the influence of its environment with an energy-preserving QND coupling. One may 
take the system to be either an oscillator (harmonic or anharmonic) or a two- level atom (or equivalently, a spin- 1/2 
system). The phase distributions associated with the quantum state for the two cases are defined. The environment is 
modeled as a bath of harmonic oscillators, and the impact of the environmental parameters is quantified for different 
initial states of the system. 

The plan of the paper is as follows. In Section II, we briefly discuss a generic energy-preserving QND system in the 
context of open systems [36j . The bath is taken to be initially in a squeezed thermal state, from which the common 
thermal bath results may be easily extracted by setting the squeezing parameters to zero. In Section III, we define the 
phase distribution for an oscillator system, following Agarwal el al. |35|. In Section IIIA in particular, we consider a 
harmonic oscillator system in QND interaction with its environment [261 ] . We study two different initial conditions, 
of the system starting (1) in a coherent state and (2) in a squeezed coherent state. In Section IIIB, we study the case 
where the system is an anharmonic oscillator, which could arise, for example, from the interaction of a single mode 
of the quantized electromagnetic field with a Kerr medium 37, 38]. This Hamiltonian can be expressed in terms of 
the generators of the group SU(1,1). Using the positive discrete series representation of this group, we construct its 
phase distribution and study it for two different situations: (1) when the system is initially in a Kerr state, and (2) 
when it is initially in a squeezed Kerr state [13]. In Section IV we consider the phase distribution for a two-level 
atom, extensively used as a model system in quantum computation [3^ . [40l l4l| . Following the phase distribution of 
angular momentum systems introduced by Agarwal and Singh [42j . we construct and study the phase distribution of 
the system for three different initial conditions of the system, starting (1) in a Wigner-Dicke state [43|, which is the 
atomic analogue of the standard Fock state [44J, (2) in an atomic coherent state, which is the atomic analogue of the 
Glauber coherent state [44| , and (3) in an atomic squeezed state [42l |4{| . In Section V we present our conclusions. 



II. GENERIC QND OPEN SYSTEMS 

We consider the following Hamiltonian describing the interaction of a system with its environment, modeled as a 
reservoir of harmonic oscillators, via a QND type of coupling [36] ]: 



H — Hs + Hr + Hsr 

l. -4- H c a,, (hu 4- h\ ) + H I 



k k 

Here Hs, Hr and Hsr stand for the Hamiltonians of the system, reservoir and system- reservoir interaction, respec- 
tively. Hs is a generic system Hamiltonian which we will specify in the subsequent sections to model different physical 
situations. b\, b k denote the creation and annihilation operators for the reservoir oscillator of frequency lju, gu stands 
for the coupling constant (assumed real) for the interaction of the oscillator field with the system. The last term on 
the right-hand side of Eq. (1) is a renormalization inducing 'counter term'. Since [Hs,Hsr] — 0, the Hamiltonian 
(1) is of QND type. The system plus reservoir composite is closed obeying a unitary evolution given by 

P {t) = e-i m p(0)ei Ht , (2) 

where 

p(0)=p s (0)^(0), (3) 
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i.e., we assume separable initial conditions. The reservoir is assumed to be initially in a squeezed thermal state, i.e., 
it is a squeezed thermal bath, with an initial density matrix Pr(0) given by 



fa(0)=S(r,*)&fc5t(r,*), 



where 



pt h = n [! 

k 



e -f3hu) k ~\ e -f3hu] k b' k b k 



(4) 
(5) 



is the density matrix of the thermal bath, and 

S(r k) <Z>k) = exp 



(b 2 P 2 
h 2 6 2 6 



(6) 



is the squeezing operator with r k , &k being the squeezing parameters |46j ]. We are interested in the reduced dynamics 
of the 'open' system of interest S, which is obtained by tracing over the bath degrees of freedom. Using Eqs. (fT]), ([3]) 
in Eq. (|2j) and tracing over the bath variables, we obtain the reduced density matrix for S, in the system eigenbasis, 
as 36] 



Pnmit) 



e -HE n -E m )t e 



x |(e Mfct - 1) cosh(r fc ) + {e~ iUki - 1) sinh(r fe )e j2 



x exp 



From ([7]) we obtain the master equation as 
where 



~{E n - E m ) + ir](t)(E 2 n - E 2 m ) - (E n - E m fj(t) 



Pnm(O)- 



PnmWi 



and 



2 fc 71 w fc V 2 y 



For the case of an Ohmic bath with spectral density 

I(u) = 2° we -<"M 

7T 

where 70 and w c are bath parameters, r](t) and 7 (i) can be evaluated [36[ and we quote the results: 

r,(t) = -^tan-^Wct), 

7T 



and 



7 (i) = cosh(2r) ln(l + w 2 t 2 ) - — sinh(2r) In 
2n 47r 

- — sinh(2r) ln(l + 4a 2 u> 2 ), 

47T 



(l + 4i^-a) 2 ) 
(l+c c 2( i -2a)2) 2 



(7) 



(8) 



(9) 



(10) 



(11) 



(12) 



(13) 
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at T = 0, with t > 2a; 



lit) 



irhiOr 



cosh(2r) 



2w c itan 1 (o; c i)+ln 



1 + LO%t 2 



■— — smh(2r) 
ZirntUr 



4w c (t — a) tan 1 (2w c (i — a)) 



-4aj c (t — 2a) tan 1 (u> c (t — 2a)) + 4aw c tan 1 (2aw c ) 



In 



[l+cjl(t-2af 
[l+4a;2(i-a) s 



In 



1 



1 + 4a 2 u 



(14) 



for high T, and again with t > 2a. Here we have taken, for simplicity, the squeezed bath parameters as 

cosh (2r(w)) = cosh(2r), sinh (2r(w)) = sinh(2r), 
$(w) = au>, 



(15) 



where a is a constant depending upon the squeezed bath. We will make use of Eqs. ©, (jTUJ) , (fl"2)) . (fl"3|) and (jT4j) in 
the subsequent analysis. Note that the results pertaining to a thermal bath can be obtained from the above equations 
by setting the squeezing parameters r and $ to zero. 

III. QUANTUM PHASE DISTRIBUTION FOR AN OSCILLATOR SYSTEM 

As discussed in the Introduction, it is more convenient to deal with the quantum phase distribution than a hermitian 
quantum phase operator. Following Agarwal et al. [351 ] we define a phase distribution V(0) for a given density operator 
p as 



27r P m > net 



(n—m)9 



(16) 



m,n— 



where the states \9) are the eigenstates of the Susskind-Glogower [29| phase operator corresponding to eigenvalues of 
unit magnitude and are defined in terms of the number states |n) as 



(17) 



n=0 



The sum in Eq. dTBJ) is assumed to converge. The phase distribution is positive definite and normalized to unity. 



A. System of a harmonic oscillator 



For the case where the system S is a harmonic oscillator with the Hamiltonian 

1 



Ha = Tilo a T a + — 



(18) 



the number states serve as an appropriate basis for the system Hamiltonian and the system energy eigenvalue in this 
basis is 



E n = fioj yn + i 
Using this in Eq. ([7]) we obtain 

p s (i) = e -i^(n-m)t e i(hu) 2 (n~m)(n+rn+l)ri(t) e -(hu) 2 (n-m) 2 ~/(t) pS ^q-j 



(19) 



(20) 
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where 77(f) and 7(f) are as in Eqs. ([9]) and (|10p . respectively. 

The Hamiltonian described here has been used by Turchette et al. [26] to describe an experimental study of the 
decoherence and decay of quantum states of a trapped atomic ion's harmonic motion interacting with an engineered 
high-temperature 'phase reservoir', which is simulated by random variations in the trap frequency - changing the 
phase of the ion oscillation without changing its energy, i.e., adiabatically modulating the trap frequency. In such 
a system it would be interesting to construct the quantum phase distribution associated with the reduced density 
matrix of the system and obtain the dynamics of the phase distribution as the system evolves under the influence of 
its environment. Equation (|20|) when substituted in Eq. (|16[) provides us with the phase distribution of the harmonic 
oscillator system interacting with its environment via a QND type of interaction. 

Now we obtain the phase distributions for some physically interesting initial conditions of our harmonic oscillator 
system S. 



1. System initially in a coherent state 



The initial density matrix of the system is 

P s (0) = \a)(a\, (21) 

where 

a = \a\e t6 ° (22) 

is a coherent state [13]. Thus the initial density matrix in the system basis is 

rfU(O) = (n\a)(a\m). (23) 

Now making use of the expansion of the coherent state in terms of the number states we get 

Ul« -id 2 . „ 

e ^ e - m9 °. (24) 

It is to be noted that each of the diagonal elements of the above density matrix (|2"3"|) is given by a Poisson distribution. 
Using Eq. (I24|) (and its complex conjugate) in Eq. ([23]) . substituting it in Eq. (120]) . and then using Eq. (|20|) in Eq. 
(|16[) we obtain the phase distribution as 

1 ^® I I TL~\~Tt\ 

■p(Q) = _ I"! p -\a\ 2 p -i(m-n)(6-8 a ) p -iu(m-n)t 

27T ^ 

m,n— U 

xe i(M 2 (m-n)(™+™+lh(t) e -(M 2 (™-"») 2 7(*) _ (25) 

Figure 1 depicts the behavior of the quantum phase distribution, V(9) given by Eq. (|25p . as a function of 9 
(in radians) as it evolves under different environmental conditions. It can be clearly seen that in comparison with 
the unitary evolution (continuous curve), as the temperature T increases, the phase distribution broadens thereby 
indicating increasing phase diffusion with T. The phase diffusion also increases with an increase in the value of the 
squeezing parameter r, defined by (|15p . as is evident from a comparison of the large-dashed and the dot-dashed curves. 
Also by comparing the dot-dashed and double dot-dashed curves indicating the same environmental conditions but 
different evolution times t, it can be seen that an increase in exposure time to the environment causes a corresponding 
increase in phase diffusion. The broadening of the curves in all cases takes place in such a fashion that the normalization 
of the phase distribution function is preserved. For all the figures in this paper, we have set u) — 1.0, ui c = 100, \a\ 2 
= 5, a = 0.0, and 9 (Eq. (22])) = 0. 

2. System initially in a squeezed coherent state 
The initial density matrix of the system is 

p s (0) = \t,a)(a,S;\, (26) 
where the squeezed coherent state is defined as [I?} 

\£,a)=S(OD(a)\0). (27) 
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FIG. 1: Quantum phase distribution, 'P(O) given by Eq. ]25H , for a harmonic oscillator initially in a coherent state, as a function of (in radians), 
for different environmental conditions and evolution times. The parameters have been taken as u) — 1.0, uj c — 100, | ct | 2 — 5, a — 0.0, 70 — 0.0025. 
The small-dashed and the large-dashed curves are for temperatures T (in units with h = ks = 1) = and 300, respectively, with an environmental 
squeezing parameter r — 2 and at an evolution time t — 0.1. The dot-dashed and the double dot-dashed curves are at evolution times t — 0.1 and 
0.2, respectively, for T — 300 and r — 1. The continuous curve represents unitary evolution (70 = 0). 



Here S denotes the standard squeezing operator and D denotes the standard displacement operator [47| . The initial 
density matrix (|26[) in the system basis is 



/Cn(0) = Hp s (0)|n} 



g 4(m-n) ( tanh ( ri ^<^+i) 



2^ 



'\fm\n\ cosh(n) 



cxp [-|a| 2 (l - tanh(r 1 )cos(26» - if)))] 



•^/sinh(2ri) 



a e 



i(fo-f) 



^/sinh(2ri) 



(28) 



where £ = r\e % ^ . Here is a Hermite polynomial. Using Eq. 

obtain the phase distribution as 



in Eq. (|20p and substituting it in Eq. (fH)|) we 



no) = 7T- 



i\ n—m 



(tanh(ri))" 



2tt ^ 

m.n— 



v /^! cosh(ri) 



x exp [-|a| 2 (l - tanh(ri) cos(26» - tp))] 



ae 



•y/sinh(2ri) 



a e 



i(»o-$) 



•^/sinh(2ri) 



^ ^—iu(m—n)t^i(huj) 2 (rn—n) (n+m+l)r/(£) ^ — (huj) 2 [n— ra) 2/ ^{t) 



(29) 



Figure 2 clearly indicates an increase in phase diffusion, corresponding to a broadening of the phase distribution 
curve, with an increase in T or bath squeezing parameter r or evolution time t, as was the case in Fig. 1. An 
interesting difference can be seen in the phase distribution curves corresponding to unitary evolution (continuous 
curves) in Figs. 1 and 2, viz., the continuous curve in Fig. 2 is more tilted than that in Fig. 1. This is due to the 
squeezing inherent in the initial state of the system (|26p which is quantified by the parameters n and ip. 
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FIG. 2: Quantum phase distribution, V{9) given by Eq. J29L for a harmonic oscillator initially in a squeezed coherent state, as a function of 9 (in 
radians), for different environmental conditions and evolution times. The parameters have been taken as u) — 1.0, H 2 = 5, a = 0.0, 7o = 0.0025: 
ri —0.5, and ip — 7r/4 (ri and ip arc the system squeezing parameters |26j). The small-dashed and the large-dashed curves are for temperatures 
T (in units with h = ks = 1) = and 300, respectively, at an environmental squeezing parameter J15t r — 2 and evolution time t — 0.1. The 
dot-dashed and the double dot-dashed curves are at evolution times t — 0.1 and 0.2, respectively, with T — 300 and r — 1. The continuous curve 
represents unitary evolution (70 — 0). 



B. System of an anharmonic oscillator 



Here we take up the case where the system is modelled as an anharmonic oscillator with the Hamiltonian 



Hs — froj 



(30) 



This has been used, for example, in studies related to a non-absorbing Kerr medium interacting with a single mode 
of the quantized electromagnetic field [13, [38[ . In such a context A is related to the third-order susceptibility of the 
Kerr medium (48|. The above Hamiltonian can be expressed (up to constant factors) in terms of the generators K , 
K + and K_ of the SU(1, 1) group. These generators have the following bosonic representation: 



K = \{<Ja + oat), K + = i( a t) 2 , K_ = 1(a) 2 . 



In terms of these generators, Eq. ([50)) can be expressed as 



H s = 2UuK + 2HXK+K. 



(31) 



(32) 



We make use of the unitary irreducible representations of the group SU (1,1) known as the positive discrete series 
V + (k), where k is the so called Bargmann inde x such that the eigenvalue of the Casimir operator of the group is 



k(k — 1) . This gives the value of k to be \ or | 50, (51( ■ The case of fc = \ marks the even sector of the representation 
with the vacuum state in this representation coinciding with the vacuum state of the harmonic oscillator, while the 



case of k 



marks the odd sector of the representation. Thus the even and the odd sectors of the representation 



together span the number states [52j. The basis for this representation obeys the following properties: 

Kolrrijk) = (m + k)\m, k), 



K+\m,k) = y/[{m + l){m + 2k)]\m + l,k), 
k) = ^[m(m + 2k - 1)] |m - 1, k) , 



(33) 



where m = 0,1,2,.... Using the above properties of the generators, the action of Hs (|30|) on the basis of this 
representation is found to be 



Hs\m,k) — 2h [uj(m + k) + \m(m + 2k — 1)] \m, k) 
= E mk \m,k). 



(34) 
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We use the above equation in Eq. (JT]) to obtain the reduced density matrix of the system S in the system basis |m, fc) 
as 

s fi\ — 2i(m — n)[u+\(m+n+2k— l)]t 

Pm k ,n h \ L ) — e 

x 4ih 2 (m-n)[u+\(m+n+2k-l)][u(n+m+2k)+\(n 2 +m 2 + (2k-l)(m+n))]r](t) 

xe -4h 2 (m~n) 2 [u+\(m+n+2k-l)] 2 ~f(t)pS ^ j-g^ ^35^) 

Now let us consider some physically interesting initial conditions for our anharmonic oscillator system S. 

1. System initially in a Kerr state 

An initial Kerr state \iPk) [131 can be obtained as a result of an interaction of the usual coherent state of the 
electromagnetic field with a nonlinear Kerr medium mediated by the Hamiltonian Hs given in Eq. (|30|) . This state 
is defined in terms of the number states as 



\M = ^2q n \n), (36) 



71 

where 

V™ _l„|2 



q n = ^=e^ e^^-V. (37) 



Here \n) represents the usual number state and \ — where A is as in Eq. (j3"0)) , L is the length of the medium and 
v is the speed of light in the Kerr medium in which the interaction has taken place. Thus the initial system density 
matrix is 

P S m k ,n k (°) = {m,k\ip K )(^K\n,k) 
= Q.2mlt n for k = 1 

3 

= 92m+l<72«+l for fc = 4' ( 38 ) 

because the state |m, fc) with fc — i represents an even number state while the state |m, fc) for fc = | represents an 
odd number state. The phase distribution is obtained by substituting Eq. (|3"5|) in Eq. (133)) and then in Eq. J16[) . 
making use of the fact that for the positive discrete series representation of the group SU(1, 1), the even and the odd 
sectors together span the number states: 

V(9) = — V g 2mg * e i2(n-m)e e -2i(m-n)[ai+A(m+n-|)]t 
m,n— 

x e <lih 2 (m-n) [ui+ X(m+n~ \ )] [u(n+m+ \ )+X(n 2 +m 2 - \ (m+n))]?)(t) 
x e -AK 2 (m-n) 2 [u+\{m+n- | )] 2 -t(t) 

i J_ \^ Oo , „* »2(n-m)fl -2i(m-n)[o)+A(m+n+i)]t 
xe 4ift 2 (m-n)[w+A(m+n+i)][w(n+m+f ) + A(n 2 +m 2 + i(m+n))]jj(i) 

xe -4R 2 (m-n) 2 [^+A(m+n+i)] 2 7 (t)^ (39) 

?2mi 92m+i can be obtained from Eq. (|3T|) . 

Figures 3 and 4 represent the evolution of the quantum phase distribution, 'P(O) given by Eq. (|39|) , as a function of 
for an anharmonic oscillator system (|30j) starting from an initial Kerr state (|36jl . While Fig. 3 represents the evolution 
for a fixed evolution time but different environmental conditions, Fig. 4 represents different evolution times under the 
same environmental conditions. From Fig. 3 it is evident that increasing the temperature T causes a broadening of 
the phase distributions. Increased phase diffusion also results from an increase in environmental squeezing r. Figure 
4 clearly shows that with an increase in the evolution time t, i.e., an increase in exposure to the environment, the 
quantum phase distribution shifts as well as diffuses. A similar conclusion was obtained by Agarwal et al. [35j for an 
analogous situation studied under unitary evolution. 
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FIG. 3: Quantum phase distribution, "P(9) given by Eq. 11391) , for an anharmonic oscillator initially in a Kerr state, as a function of (in radians), 
for different environmental conditions at a fixed time of evolution. The parameters have been taken as 70 — 0.0025, | ct | 2 — 5, to — 1.0, x — ^ — 
0.02 and evolution time t — 0.1. The small-dashed and the large-dashed curves are for the bath squeezing parameter r — and 2, respectively, at 
T (in units with h = ks = 1) = 50. The dotted curve is for T — and r — 2. The continuous curve represents unitary evolution (70 — 0). 

P (6) 
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FIG. 4: Time evolution of the quantum phase distribution, "P(#) given by Eq. ]39t , for an anharmonic oscillator initially in a Kerr state, as a 
function of 8 (in radians), for different evolution times under fixed environmental conditions. The parameters have been taken as 70 — 0.0025, \a\ 2 
— 5, uj — 1.0, x — ^ — 0.02 (as in Fig. 3), T — 0, and r — 2. The dotted curve is for an evolution time t — 0.1, the small-dashed curve is for t — 
0.5, and the large-dashed curve is for t — 1.0. 



2. System initially in a squeezed Kerr state 
A squeezed Kerr state is defined in terms of the number states as 

\ipSK) = ^2s m \m). 



(40) 



Thus the initial system density matrix in the system basis |m, k) (134|) is 

Pm fc ,™ fe (0) = {m,k\^ SK ){%p SK \n,k) 
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S2rnS* 2n for k = \ 



S2 m +l4n+l for k = J- 



Here 



S2m — <l2pG2m2p{z), 



(41) 
(42) 



and 



S2m+1 — g2p+l<J2m+12p+l(^); 



(43) 



with z = r^e 4 *, and G mp (z) 



G 



(m\S(z)\p), where S(z) is the usual squeezing operator, is given by [5G 
(-I)p ((2p)\(2m)V 



2m2p 



p\m\ 



cosh(ri) 



exp (i(m — p)(j)) 



( tanh(ri] 



(m+p) 



Ft 



Similarly, G 2 m+i2p+i(^) is given by 



G 



2m+12p+l 



{-l) p f(2p+ l)!(2m+ 1)! 



p\m\ 



cosh (ri ) 
/tanh(n)\ (m+p) 



J 



'2' (sinh(n)) 2 

exp (i(m — p)</>) 
3 1 



(44) 



(sinh(ri)) 2 



(45) 



Here Ff is the Gauss hypergeometric functio n 15411 . We make use of even p in Eqs. (|42| . (|44]) , and odd p in Eqs. 



(|43| . (|45|) . because as has been pointed out in [53j, G mn is nonzero only for m, n either both even or both odd. Since 
m is even in (|4"2"|) . it follows that p too should be even, and similarity for (|4"3")l where m is odd. For convenience, it is 
sometimes assumed that </> is zero and z — r\ is real. Using Eqs. ([42]) . (|43|) . (|44|) and (|45]) in Eq. (|4T|l . substituting it 
in Eq. (|35p and then in Eq. (|16p we obtain the phase distribution function as 



V(9) = 



1 „ „* J2(n-m)6 -2i(m-n) L+A(m+ra-|)lt 

Z7T * — ' 

m,n— 

x g4i^ 2 (m-n) [a>+ \(m+n- \ )] [o)(n+m+ i )+A(n 2 +m 2 - \ (m+n))]?)(t) 
xe -4fi 2 (m-n) 2 [^+A(m+n-I)] 2 7 (t) 



,J_ * i2(n-m)9 -2i(m-n)[u>+A(m+n+|)]t 

+ 27T * 2m + l6 2n+l e e 

m,n— 

xe 4i^ 2 (m-n)[w+A(m+n+i)][w(n+m+f ) + A(n 2 +m 2 + i(m+n))]jj(t) 
xe -4fi 2 (m-n) 2 [w+A(m+n+i)] 2 7(t) 



(46) 



Here S2m is as in Eq. (|4"2"| and S2m+i is as in Eq. (|4"3"1) . 

Figure 5 depicts the evolution of the quantum phase distribution, V{&) given by Eq. ([46]) . as a function of 9 for 
an anharmonic oscillator system (|30[) starting from an initial squeezed Kerr state (|4U|) . The environmental effects 
are clearly depicted in that an increase in temperature T and squeezing parameter r results in the broadening of the 
phase distribution indicating increased phase diffusion. In the same figure, a drastic influence of the environment on 
the unitary behavior can be seen from the sharp fall in the amplitude of the phase distribution with the inclusion of 
environmental effects. A comparison between the unitary evolutions (continuous curves) of Figs. 3 and 5 highlights 
the difference in the initial conditions of the system depicted in these curves. The peak amplitude of the unitary 
evolution (continuous curve) is greater in case of a system initially in a squeezed Kerr state (Fig. 5) than that in a 
Kerr state (Fig. 3). This is indicative of the additional squeezing in the initial state for Fig. 5. The corresponding 
narrowing of the peaks in Fig. 5 is due to the fact that the phase distributions are normalized. Note that the multiple 
peaks are a common feature of the quantum phase distributions of the anharmonic oscillator system (|30[) as opposed 
to the single peaks of the harmonic oscillator system (jT5J) . 
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FIG. 5: Quantum phase distribution, 'P(O) given by Eq. J46| l, for an anharmonic oscillator initially in a squeezed Kerr state, as a function of 9 
(in radians), for different environmental conditions. The parameters have been taken as t — 0.1, X — ^ — 0.02, 70 — 0.025, Ti — 0.4 and — (ri, 
(j> are the system squeezing parameters). The dot-dashed and the dotted curves are for T (in units with h = ks = 1) = and 100, respectively, 
with the environmental squeezing parameter r — 1. The large-dashed curve is at T — 100 and r — 0. The continuous curve represents the unitary 
evolution (70 — 0). 



IV. QUANTUM PHASE DISTRIBUTION OF A TWO-LEVEL ATOMIC SYSTEM 

In this section we discuss the case where our system S is a two-level atom with a representation of the group SU{2). 
The system Hamiltonian, to be substituted in Eq. (JTJ) , is 



huj 



(47) 



where az is the usual Pauli matr ix ( as has been used, for example, in the quantum computation models in [39l. |40L l41|| ) . 
The Wigner-Dicke states [43|, |44| [55| \j, m), which are the simultaneous eigenstates of the angular momentum operators 
J 2 and Jz, serve as the basis states for Hg and we have 



Hs\j,m) = twjm\j,m) 
= Ej. m \j,m) 



(48) 



Here —j < m < j. Using this basis and the above equation in Eq. ([7]) we obtain the reduced density matrix of the 
system as 



s ( f \ _ -iui(m-n)t i(huj) (m -n )rj(t) 
l J jm,jn\ L ) ~ e e 

,(0). 



-(ftw)^(m-n) 2 7(i) s 

x e p ifi 



(49) 



Following Agarwal and Singh [42| we introduce the phase distribution V (</>), <p being related to the phase of the dipole 
moment of the system, as 



4tt jo 

r 2TT 



(50) 



where V{4>) > and is normalized to unity, i.e., J Q v d(f>V{4>) = 1. Here Q(8, (f>) is defined as 

Q{6,4>) = (8,<p\p s \O,0), (51) 
where \6,<j>) are the atomic coherent states [IB, H3] given by an expansion over the Wigner-Dicke states [HJ as 

3 



E 

m=-j 



2j 
j + m 



(sin(e/2)) j+m (cos(d/2)y- m \j,m)e-^ j+m ^. 



(52) 
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Using Eq. (|5T]) in Eq. ([50|) , with insertions of partitions of unity in terms of the Wigner-Dicke states, we can write 
the phase distribution function as 



m - 2} + 1 r 



47T 



/ dO sin 6» V (9,(j>\j, 
Jo „ , „• 



n,m— — j 

x(j,n\p s (t)\j,m)(j,m\e,4>). (53) 

We make use of Eq. (|4"9"f and 

i 

(j, m|0, 0) = L 2 (sin(0/2))-'"+ m (cos(0/2))^™ e - 4 (-'"+ m )* (54) 



and its conjugate in Eq. (|53|) to obtain the required phase distribution for specific initial conditions of the system 5*. 
Let us now consider some physically interesting initial conditions for the two-level system S. 

A. System initially in a Wigner-Dicke state 

A Wigner-Dicke state is the atomic analogue of the Fock state [Hj]. The initial density matrix of the system S in 
this case is 

p s (0) = \j,m)(j,fh\, (55) 

which gives 

(j, n\p S (t)\j, m) = 6 n ,rn6m,m- (56) 

Using this, the phase distribution becomes 

Vtt) = ( „• ) B U +fh+l,j-m+l}. (57) 



2ir \j+m 

Here B stands for the Beta function. It is evident from Eq. ([57)1 that the phase distribution for the atomic system 
starting in a Wigner-Dicke state is uniform and is independent of any bath dynamics. Here, since we have only one 
two-level system in Eq. (|47p. j = \ and V(4>) can be seen to go over to i.e., a uniform distribution. 

B. System initially in an atomic coherent state 

An atomic coherent state is the atomic analogue of the Glauber coherent state [44| . The initial density matrix of 
the system S in this case is 

P s (0) = \a,fl(a,/3\ (58) 
yielding the following matrix element in the \j, m) basis: 

U,n\p s (t)\j,m) = 



— iuj(n — m)t i(huj) 2 (n 2 —m 2 )r)(t) 



Xe -(M (n-m) 7(*)^„| 0>) g)( a)/ g| J - sTO ). (59) 
Using Eqs. ((SH), in Eq. (p)3")) we obtain the phase distribution as 

V + 1 r M ^a ST ( 2i \ ( 2j 



j -\- n j u + m 



n,m——j 

2j — n—rn —i(n—m){3 



x (sin(6l/2)) 2j+ ™ +m (cos(6l/2)) 
x(sin(a/2)) 2 ''' +n+ " i (cos(a/2)) 23 -"- m 

— iuj(n — m)t i(huj) 2 (n 2 — m 2 )r](t) 
xe -( ?iw ) 2 (™-« l ) 2 7(t) e i(™-™)0 > (gO) 
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FIG. 6: Quantum phase distribution, "P{(p) given by Eq. H62i , for a two-level atom initially in an atomic coherent state, as a function of <ft ( m 
radians), for different environmental conditions and evolution times. The parameters have been taken as a — (3 — 7r/4 [Eq. ]58H , and 70 — 
0.025. The continuous curve and the dotted curve are for the bath squeezing parameter r — and 2, respectively, at a temperature T — and an 
evolution time t — 0.1. The small-dashed and the large-dashed curves correspond to evolution times t — 0.1 and 0.02, respectively, at T (in units 
with h = k B = 1) = 300 and r = 2. 



In the above equation, the 6 integral can be carried out to yield 

p(0) = 2j + 1 j2 C 2j 1 ^ 2j \ r ( j + ^ n+m ") + v r ( j -^ n+m ^ + 1) 



2tt ^ \j + nj\j + ml F(2j + 2) 

n,m——j 

xe-^ n - m ^(sin(a/2))^+ n + m (cos(a/2)) 2j - n - m 

^ g—iLu(n—m)t^i(hu) 2 (n 2 —m 2 )rj(t) 

Here T is the standard Gamma function. Since H$ given by Eq. (|47|) represents a single two- level atom, j = \- Eq. 
(I6ip is thus considerably simplified and we obtain the phase distribution as 



1 



2tt 



7T 



V{4>) = ^- l+'-smacos{f3 + ojt-(t))e' {hul) l{t) . (62) 



4 



It can be easily checked that this V (<f>) is normalized to unity. As can be seen from Eq. (|62|) . only j(t) plays a role in 
carrying the effect of the environment on the phase distribution. For a generic QND open quantum system described 
by (|I1), it can be shown that j(t) is the decoherence causing term (36|. Thus Eq. (|62)) is a simple and neat formula 
clearly illustrating the effect of the environment on phase diffusion. By making use of y(t) from Eqs. (|13| and (fT4| for 
T = and for high T, respectively, we find that the second term on the right-hand side of Eq. (|6"2"]) has a power-law 
decay at zero T and an exponential decay at high T, and eventually the phase distribution tends to the uniform limit 
of Thus the effect of the environment stays for a longer time at zero T as compared to that at high T. 

Figure 6 depicts the evolution of the quantum phase distribution, V{<j>) given by Eq. (|62|) . as a function of 4> ( m 
radians) for different environmental conditions. It is clearly seen that increasing the temperature T, the bath squeezing 
parameter r, and the environment exposure time t cause a broadening of the phase distribution curve, indicating an 
increase of phase diffusion. The broadening of the curves preserves the normalization of the phase distribution. 

C. System initially in an atomic squeezed state 

An atomic squeezed state [IE], [HB, HH is expressed in terms of the Wigner-Dicke states as 

10 P) = A p exp(Ojz)exp(-i^J Y )\j 1 p), (63) 
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where 

e 2e = tanh(2|C|), (64) 
with £ indicating the initial squeezing of the system. The initial density matrix of the system S in this case is 



p s (0) = \(,p)({,p\- 

Using Eq. (|53|) along with the expressions 

{j,n\p s (t)\j,m) = e-Mn-^t^fi^-^Mt) 

xe -W(n-m ) 2 T(*)( J> | C ,p)( C , p | i , m ) ) 

and 

(j,n\C,p)=A p e ne di p (^), 
[56| . where <P np (^) is the standard Wigner symbol for the rotation operator (59j : 



and 



i^i 2 = (E 

we obtain the phase distribution function as 

2 j + 1 , 



(-l) r (2j-r)!(cosh6) 
H(j +p- r)!(j -p - r)\ 



(65) 

(66) 
(67) 

(68) 
(69) 



4tt 



A p | 2 / d0 sin ^ 



2j 



2j 
j + m 



x(sin(6»/2)) 2j+ " +m (cos(6'/2)) 2j - n - m e 4( "- m)0 



In Eq. (|70|) the integral can be carried out to yield 



(70) 



V{cj>) 



2j + l 
2tt 



- 1 / 



n,m——j 



2.7 



j + n J \j + m 



2j \ * T(j + |(n + m) + l)r(j - |(n + to) + 1) 



r(2j + 2) 



xe 

xe -(r lW ) 2 (n-m) 2 7(t) e (™+m)e 



(2) (2) 



(71) 



As discussed in subsection IVB above, for a single two-level system, j = i. We take up two cases for the two values 
of p appearing in Eq. ((63]) : p — — ^ called the south pole of the phase space of the two-level system, and p—\ called 



(72) 



the north pole of the phase space [4 

a. South pole (p = — |): 
The phase distribution in Eq. (fTTj) is considerably simplified to give 

1 



b. North pole (p = i): 
Eq. (|7ip is simplified to give 



V{4>) = [l + , \^ cos(0 - wtJe-^'TW 
v y 27r L 4 cosh 8 



(73) 
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FIG. 7: Quantum phase distribution, V{4>) given by Eq. J72H , at the south pole of the phase space for a two-level atom initially in an atomic 
squeezed state, as a function of <f> (in radians), for different environmental conditions. The system squeezing parameter ]64t is taken as — - 
0.5494, the bath squeezing parameter r — 1, and 70 — 0.025. The continuous curve corresponds to a temperature T — and an evolution time t 
— 0.1, while the small-dashed and large-dashed curves correspond to evolution times t — 0.1 and 0.05, respectively, at T (in units with Ti = ks = 
1) = 300. 
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FIG. 8: Quantum phase distribution, "P{4>) given by Eq. J73}, at the north pole of the phase space for a two-level atom initially in an atomic 
squeezed state, as a function of (in radians), for different environmental conditions. The system squeezing parameter J64t has been taken as 
— - 0.5494, the bath squeezing parameter r — 1, and 70 — 0.025. The continuous curve corresponds to a temperature T — and an evolution time 
t — 0.1, while the small-dashed and large-dashed curves correspond to evolution times t — 0.1 and 0.05, respectively, at T (in units with h = \zb = 
1) = 300. 

A p is defined by Eqs. (|63[) and (|67p . and is usually fixed by normalization as in the above equations, where 
\A p \ 2 = (coshO) -1 which is equal to 1 for 8 = [56[, implying an infinite initial squeezing £ of the system (|64| . 
The expression in Eq. (|72p for the south pole in the phase space and that in Eq. (|73|) for the north pole in the 
phase space differ from each other by a sign in the second part of the expressions. The contrast between them can 
be seen clearly from Figs. 7 and 8. As seen from the figures, with the increase in temperature T or the exposure 
time to the environment t, the phase distribution curves flatten out indicating increased phase diffusion. These curves 
bring out another notable feature, viz., with the increase in bath exposure time t or temperature T, the effect of 
squeezing, indicated by the parameter r, is washed out. This behavior is analogous to the effect of squeezing in 
oscillator systems, where with the increase in bath exposure time t and temperature T, the non-stationary effects 
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introduced by the squeezed bath are washed out [60(. 

Eqs. (|72p and ([75)) are easily seen to be normalized to unity. As in subsection IVB above, the effect of the 
environment shows up in the above equations only in the function 7(f), responsible for decoherence. As was the case 
in the previous subsection IVB, from the forms of the function j(t), for an Ohmic bath, given by Eqs. (flU)) and (jT5J) 
for T = and for high T, respectively, we find that the second term on the right-hand side of Eqs. ([72]) and (|T3[) has 
a power-law decay at zero T and an exponential decay at high T, and eventually the phase distribution tends to the 
uniform limit of ^ seen in the case in IVA. Thus the effect of the environment stays for a longer time at zero T as 
compared to high T, and eventually the distribution tends to the same uniform value irrespective of the initial state 
being a coherent or a squeezed state. As pointed out in [H,[45|, H2I, the state \C,p) in |63|) has an inherent squeezing 
which is represented by £ in (|64[) and the environmental squeezing r is encapsulated in the function "f(t) given by 
(fiTj]) . Thus the results (fT2"]) and (|T3j) bring out the relative importance of the different squeezing sources and hence 
are applicable in the context of the experiment of Kuzmich et al. (23j concerning the role of the environment on the 
atomic quantum nondemolition measurements and squeezing. 

V. CONCLUSIONS 

In this paper we have analyzed the quantum phase distribution of a number of physically interesting systems 
interacting with their environment via a QND type of coupling. We have taken our system to be cither an oscillator 
(harmonic or anharmonic) or a two-level atom (or equivalently, a spin-1/2 system), and modeled the environment as 
a bath of harmonic oscillators, initially in a general squeezed thermal state, from which the common thermal bath 
results may be easily extracted by setting the squeezing parameters to zero. We have explicitly evaluated the phase 
distribution and worked out the effects of different environmental parameters on the dynamics of the system starting 
with various initial states. 

In particular, for a harmonic oscillator system in QND interaction with its environment (Section IIIA), we have 
considered two different initial conditions of the system, starting (1) in a coherent state, and (2) in a squeezed coherent 
state. The phase distribution corresponding to the unitary evolution in the second case is more tilted than that in 
the first case, which is a signature of the squeezing inherent in the initial state of the system. We have next taken 
an anharmonic oscillator (Section IIIB), which could arise, for example, from the interaction of a single mode of the 
quantized electromagnetic field with a Kerr medium, and constructed its phase distribution, again for different initial 
conditions: (1) the system starting in a Kerr state, and (2) the system starting in a squeezed Kerr state. With an 
increase in the evolution time t, indicating an increase in exposure to the environment, the quantum phase distribution 
in the first case shifts as well as diffuses. 

We have then studied the phase distribution for a discrete two- level atom (Section IV) , for different initial conditions 
of the system, starting (1) in a Wigner-Dicke state, which is the atomic analogue of the standard Fock state, (2) in 
an atomic coherent state, which is the atomic analogue of the Glauber coherent state, and (3) in an atomic squeezed 
state. In the first test case, the phase distribution is uniform and is independent of any bath dynamics. In the other 
two cases, it is seen that the effect of the environment stays for a longer time at zero temperature than that at high 
temperature, and eventually the distribution tends to the uniform value of the first case, irrespective of the initial 
state of the system. 

In all the cases considered, a broadening of the phase distribution curve, indicating an increase in phase diffusion, 
results with an increase in the bath temperature T or bath squeezing parameter r or evolution time t. The broadening 
of the curves, of course, preserves the normalization of the phase distribution. Even though each system considered 
is an 'open' system, we could make use of the underlying group symmetries of the system Hamiltonians, because of 
the QND nature of the system-environment coupling. Our quantitative results are of potential use in the analysis of 
a broad class of relevant experimental situations dealing with quantum nondemolition measurements and squeezing. 
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